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Local field potentials (LFPs) are routinely measured experimentally in brain tissue, and exhibit 
strong low-pass frequency filtering properties, with high frequencies (such as action potentials) 
being visible only at very short distances («10 pm) from the recording electrode. Understanding 
this filtering is crucial to relate LFP signals with neuronal activity, but not much is known about 
the exact mechanisms underlying this low-pass filtering. In this paper, we investigate a possible 
biophysical mechanism for the low-pass filtering properties of LFPs. We investigate the propagation 
of electric fields and its frequency dependence close to the current source, i.e. at length scales in 
the order of average interneuronal distance. We take into account the presence of a high density 
of cellular membranes around current sources, such as glial cells. By considering them as passive 
cells, we show that under the infiuence of the electric source field, they respond by polarisation, 
i.e., creation of an induced field. Because of the finite velocity of ionic charge movement, this 
polarization will not be instantaneous. Consequently, the induced electric field will be frequency- 
dependent, and much reduced for high frequencies. Our model establishes that with respect to 
frequency attenuation properties, this situation is analogous to an equivalent RC-circuit, or better 
a system of coupled RC-circuits. We present a number of numerical simulations of induced electric 
field for biologically realistic values of parameters, and show this frequency filtering effect as well 
as the attenuation of extracellular potentials with distance. We suggest that induced electric fields 
in passive cells surrounding neurons is the physical origin of frequency filtering properties of LFPs. 
Experimentally-testable predictions are provided allowing to verify the validity of this model. 

PACS numbers: 87.19. La, 87.17.Aa 



I. INTRODUCTION 

Electric fields of the brain which are experimentally observable either on the surface of the scalp, or 
by microelectrodes, are due to currents in dendrites and the soma of cortical pyramidal cells P, 01- Ex- 
perimental measurements of electric fields created in the brain distinguish three scenarios: (i) Local field 
potentials (LFPs) denote the electric potential recorded locally in the immediate neighborhood of neurons 
using microelectrodes of size comparable to the cell body, (ii) The electrocorticogram (ECoG) refers to 
measurements of the field using electrodes of a diameter of the size of about 1 mm, placed on the cortical 
surface, (iii) The electroencephalogram (EEG) is measured at the surface of the scalp using electrodes of a 
centimeter scale. In the latter case, the electric potential is recorded after conduction through ccrobrospinal 
fiuid, cranium and scalp, and corresponds to the situation where the source of the electric signals in the 
cortex is located far from the site of detection on the scalp (at a scale A « 10'^(i„„, where d„„ = 0.027 mm 
is the average distance between cortical neurons). 

By contrast with the intracellular or membrane potential, which biophysical properties have been exten- 
sively studied [^^|^, the mechanisms underlying the genesis of LFPs arc still unclear. LFP recordings 
routinely show strong high-frequency attenuation properties, because action potentials are only visible for 
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a few neurons immediately adjacent to the electrode, while low-frequency components results from large 
populations of neurons in the local network. Using LFPs, it has been shown that low-frequency oscillations 
(0-4 Hz) have a large scale coherence, while the coherence of higher frequency (20-60 Hz) oscillations is 
short-ranged. In anesthetized animals, it was shown that oscillations of up to 4 Hz have a coherence range 
in the order of several millimeters, while oscillations of 20-60 Hz have a coherence of sub-millimeter range 
The local coherence of high frequency oscillations was also shown in the visual system of anesthetized 
cats, where gamma-oscillations (30-50 Hz) appear only within restricted cortical areas and time windows 

I3 • The same difference of coherence between low and high frequency oscillations was also demonstrated 
in non-anesthetized animals, respectively during sleep in a state of low- frequency wave and waking |10| . 
Similar findings have been reported for human EEG [ij] . 

A full understanding of the mechanisms underlying the genesis of EEG and LFP signals is required to 
relate these signals with neuronal activity. Several models of EEG or LFP activity have been proposed 
previously (e.g., see [11 [13 13 llS 13 113 ) ■ These models always considered current sources embedded 
in a homogeneous extracellular fluid. In such homogeneous media, however, there cannot be any frequency 
filtering property. Extracellular space consists of a complex folding of intermixed layers of fluids and 
membranes, while the extracellular fluid represents only a few percent of the available space 0,1131 . Due 
to the complex nature of this medium, it is very difficult to draw theories or to model LFPs properly, and 
one needs to make approximations. In a previous paper [l^ . we considered current sources with various 
continuous profiles of conductivity according to a spherical symmetry, and showed that it can lead to 
low-pass frequency filtering. This showed that strong inhomogeneities in conductivity and permittivity in 
extracellular space can lead to low-pass frequency filtering, but this approach was not satisfactory because 
high-pass filters could also be obtained, in contradiction to experiments. In addition, this model predicted 
frequency attenuation which was not quantitative, as action potentials were still visible at 1 mm distance, 
which is in contrast to what is observed experimentally. 

In the present paper, we go one step further and consider an explicit structure of extracellular space, in 
which we study the interaction of the electric field with the membranes surrounding neurons. Neurons are 
sourrounded by densely packed membranes of other neurons and glial cells (19f2Q..2L] . a situation which we 
approximate here by considering a series of passive spherical membranes around the source, all embedded 
in a conducting fluid. We show that low-pass frequency filtering can be determined by the membranes of 
such passive cells and the phenomenon of electric polarisation. 



II. GENERAL THEORY 



In this section, we describe the theoretical framework of the model. We start by outlining the model 
(Section |II A\ ). where a simplifled structure of extracellular space is considered. We also describe its 
main simpliflcations and assumptions. Next, we discuss the physical implementation of this model (Sec- 
tion jll Bp . while the propagation of the electric fields will be analyzed more formally in the next sections 
(Sections HEOIIIE])- 



A. Model of extracellular space 



In Ref.jlg, we considered a model where the electrical properties of extracellular space are described 
by two parameters only: conductivity a and permittivity e. Both, a and e were considered to vary with 
location according to some ad hoc assumptions. Kj, the frequency uj component of the electric potential V 
obtained by a Fourier transformation of the potential as function of time, was shown to obey the equation 



(1) 



3 



The physical behavior of the solution is essentially determined by the expression 1 + iu-. When a strong 
inhomogeneity occurs, one may have >> 1. This means that a strong phase difference may exist between 
electric current and potential, i.e. a large impedance. In a neurophysiological context, such behavior is 
likely to occur at the interface between the extracellular fluid (high conductivity) and the membranes of cells 
(low conductivity). This situation was considered previously in the context of a simplified representation of 
extracellular space, in which the inhomogeneities of conductivity were assumed to be of spherical symmetry 
around the source Moreover, the equation above assumes that the charge density of the extracellular 
medium is zero when its potential is zero, which is not entirely true since there is an excess of positive 
charges on the exterior surface of neuronal membranes at rest. 

In the present work, we consider an explicit structure of extracellular space to more realistically account 
for these inhomogeneities of conductivity. The extracellular space is assumed to be composed of active cells 
producing current sources (neurons), and passive cells (glia), all embedded in a conducting fluid. Neurons 
are characterized by various voltage-dependent and synaptic ion channels, and they will be considered 
here as the sole source of the electric field in extracellular space. On the other hand, glial cells are very 
densely packed in interneuronal space, sometimes surrounding the soma or the dendrites of neurons |20ll21| . 
Glial cells normally do not have dominant voltage-dependent channel activity, and they rather play a role 
in maintaining extracellular ionic concentrations. Like neurons, they have an excess of negative charges 
inside the cell, which is responsible for a negative resting potential (for most central neurons, this resting 
membrane potential is around -60 to -80 mV). They will be considered here as "passive" and representative 
of all non-neuronal cell types characterized by a resting membrane potential. We will show that such passive 
cells can be polarized by the electric field produced by neurons. This polarization has an inertia and a 
characteristic relaxation time which may have important consequences to the properties of propagation of 
local field potentials. These different cell types are separated by extracellular fiuid, which plays the role 
of a conducting medium, i.e. allows for the flow of electric currents. In the remainder of this text, we will 
use the term "passive cell" to represent the various cell types around neurons, but bearing in mind that 
they may represent other neurons as well. 

Another simplification is that we will consider these passive cells as of elementary shape (spherical 
or cubic). Under such a simplification, it will be possible to treat the propagation of field potentials 
analytically and design simulations using standard numeric tools. Our primary objective here is to explore 
one essential physical principle underlying the frequency-filtering properties of extracellular space, based 
on the polarization of passive membranes surrounding neuronal sources. We assume that such a principle 
will be valid regardless of the morphological complexity and spatial arrangement of neurons and other cell 
types in extracellular space. As a consequence of these simplifications, the present work does not attempt 
to provide a quantitative description but rather an exploration of first principles that could be applied in 
later work to the actual complexity of biological tissue. 

The arrangement of charges in our model is schematized in Fig. ^^J, where we delimited 5 regions. The 
membrane of the passive cell (region 3) separates the intracellular fluid (region 5) from the extracellular fluid 
(region 1), both of which are electrically neutral. The negative charges in excess in the intracellular medium 
agglutinate in the region immediately adjacent to the membrane (region 4), while the analogous region at 
the exterior surface of the membrane (region 2) contains the positive ions in excess in the extracellular 
space. This arrangement results in a charge distribution (schematized in Fig. 03]) which creates a strong 
electric field inside the membrane and a membrane potential. 

The behavior of such a system depends on the values of conductivity and permittivity in these different 
regions (they are considered constant within each region). The extracellular fiuid (region 1) has good 
electric conductance properties. We have taken as conductivity tJi = 4 Ohm~^meter~^, consistent with 
biological data a = 3.3 — 5 Ohm~^meter~^, taken from measurements of specific impedance of rabbit 
cerebral cortex This value is comparable to the conductivity of salt water {agw = 2.5 Ohm^^ meter"^). 
The permittivity is given by ei =70 eq, corresponding to salt water. Here eo = 8.854 x 10~^^ Farad/meter 
denotes the permittivity of the vacuum. In region 2, to the our best knowledge, there are no experimental 
data on conductivity close to the membrane. We have chosen the values of a2 ~ 0.7 x 10~^ Ohm^^metcr"^ 
and 62 = 1.1 X Farad metcr"^ k, 12 eo for region 2. Such a choice is not inconsistent with biological 
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FIG. 1: Scheme of charge distribution around the membrane of a passive ceU. A. Charge distribution at rest. The 
following regions are defined: the extracellular fluid (region f), the region immediately adjacent to the exterior of 
the membrane where positive charges are concentrated (region 2), the membrane (region 3, in gray), the region 
immediately adjacent to the interior of the membrane where negative charges are concentrated (region 4), and 
the intracellular (cytoplasmic) fluid (region 5). B. Schematic representation of the charge density as a function of 
distance (along the horizontal dotted line in A). C. Redistribution of charges in the presence of an electric field. 
The ions move away from or towards the source, according to their charge, resulting in a polarization of the cell. 
D. Schematic representation of the charge density predicted from C. 



observations. First, electron microscopic photographs taken from the region near the membrane reflect very 
httle hght, which hints to quite low conductivity compared to the conductivity of region 1. We consider 
it as plausible that permittivity in region 2 should be smaller than in region 1. Our choice of a2 and 62 
corresponds to a Maxwell time Tm yielding a cut-off frequency fc ~ 100 Hz, which was also the choice 
given in a previous study investigating composite materials |23j . 

For passive cells, we neglect ion channels and pumps located in the membrane, which is equivalent to 
assume the absence of any electric current across the membrane. Therefore, Region 3 has zero conductivity 
perpendicular to the membrane surface. The capacity of a cellular membrane has been measured and is 
about C = 10^^ Farad/meter^ Approximating the membrane by a parallel plates capacitor (with 
surface S and distance d, obeying C = tS/d), one estimates the electric permittivity of membrane to 
£3 = 10~^" Farad/meter. Hence we used the parameters CT3 = 0, and £3 = 12 eo- 

Thus, the basic idea behind the model is as follows. As represented in Fig. PP, we consider a single 
spherical passive cell under the influence of an electric field. The electric field will induce a polarisation 
of the cell by reorganizing its charge distribution (Fig. Qp-D]). This polarisation will create a secondary 
electric field, with field lines connecting those opposite charges. It is customary notation to call the original 
electric field the source field, or the primary field, while the field due to polarisation is called the induced 
field, or the secondary field. The physical electric field is the sum (in the sense of vectors) of both, the source 
and the induced field. This induced field will be highly dependent on frequency, for high frequencies, the 
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"inertia" of charge movement in regions 2 and 4 will limit such a polarisation, and will reduce the effect of 
the induced field. This phenomenon is the basis of the model of frequency-dependent local field potentials 
presented in this article. 

B. Physical implementation of the model 

We introduce here the toolbox of physics used for the construction of our model, in which we consider 
neurons as the principal sources of electric current. In electrodynamics, one distinguishes between a current 
source and a potential source. A current source, e.g. in an electric circuit, means that if the electric current 
is used to do some work, the source maintains its level of current. Likewise, a potential source maintains 
its electric potential. Although neurons are generally considered as current sources, here we have chosen to 
consider them as potential sources, for the following reasons. First, these two types of sources are equivalent 
for calculating extracellular potentials. The electric potential within a given region D depends on the limit 
conditions at the border of this region, and not on the type of source (current vs. potential source). Second, 
it is much easier to calculate the potential using limit conditions on the potential (Dirichlet conditions) 
than using limit conditions with currents (Neumann conditions). Third, potentials are better constrained 
from experiments and their range of values is better known than currents (for example the amplitude of 
membrane potential variations in neurons is of the order of 10-20 mV for subthreshold activity, and of 
about 100 mV during spikes). 

Let us suppose that we are given a source of electric field, representing a neuron with some open ion 
channels. The motion of ions through those channels gives a current density, and also creates a distribution 
of charges (ions). According to the Maxwell equations, a given distribution of currents and charges (plus 
information on the polarisation and induction properties of the medium) determines the magnetic field, 
the electric field and hence the electric potential. In a first step, we consider the effect of the electric field 
on a single passive cell. 

The source electric field is the origin of two physical phenomena. First, the electric field exerts a force on 
charged particles (ions) and thus creates a motion of charge carriers, i.e. an electric current. Second, the 
electric field creates a displacement of charges in the borderline region of the membrane of the passive cell 
(region 2). This displacement creates an induced potential due to polarisation. However, this polarisation 
is not instantaneous, due to the "inertia" of charge movements. The charges on the membrane move 
relatively slowly, which will be responsible for a slow time dependence of the polarisation, before reaching 
equilibrium. The characteristic time scale of such a process is given by Maxwell's relaxation time Tm 
(see Eq. [7| below), which depends on the properties of the medium, like resistivity and permittivity. The 
temporal behavior of the source is dynamic, which can characterized by some characteristic time scale Tg. 
E.g., during the creation of an action potential, ion channels open and ionic currents flow, the electric 
field and potential changes. After a certain time the neuron source goes back to its state of rest. This 
means that the electric current will vanish after some delay. For example, in a typical neuron, the width 
of the action potential is typically of the order of 2 msec. The polarisation and relaxation dynamics of the 
passive cell depends on both time scales, Ts and Tm- This situation holds in the case of an ideal dielectric 
medium (conductivity zero, for example pure water). In the extracellular medium, however, conductivity is 
non-zero. As a consequence, the polarisation potential and the current distribution will mutually influence 
each other. 

Thus, we face the question: How do we quantitatively calculate the electric field and potential in such a 
medium? And what is its time-dependence and frequency-dependence? Looking at the temporal behavior, 
there are two different regimes from the biological point of view. First, there is the transitory regime which 
describes the short period of opening and shutting down of the source and its response of the medium with 
some delay. Second, there is the so-called asymptotic or permanent regime, where no current flows. These 
regimes will be considered in the next sections. 
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C. Asymptotic behavior in region 1 

Let us start by considering the behavior in region 1 when the system has settled into an asymptotic 
regime. It means that charges move and after some time attach to the surface of the passive cell membrane 
(region 2). We call this a stationary or equilibrium regime. In this region, the electric properties are 
characterized by permittivity ei and conductivity CTi . We assume that these electric properties are identical 
everywhere inside region 1, i.e. permittivity ei and conductivity ai are constant. 

To describe this situation, we start by recalling the fundamental equations of electrodynamics. Gauss' 
law, which relates the electric field E to the charge density p, reads 

V • {e,E) =p. (2) 

Moreover, because cri is constant, there is Ohm's law, which relates the electric field E to the current 
density j, 

J^^iE. (3) 
From Maxwell's equations one obtains the continuity equation 

V.J.I^O^ (4, 

Using Eqs. j2l4| and recalling that ei and cti are both constant in region 1 implies the following differential 
equation for the charge density, 

dp <Ji . . 

A particular solution to this equation requires to specify boundary conditions and initial conditions. Here 
the boundary conditions arc such that the electric field (the component of the field perpendicular to the 
surface) on the surface of the source and on the surface of the passive cell is given (Dirichlct boundary 
conditions). The general solution of Eq. is 

p{x, t) = p{x, 0) exp[-—t] = p{x, 0) exp[-i/Tfl/J , 
£1 

p{x, t) — H^oc , (6) 

i.e., with increasing time the charge distribution goes exponentially to zero. The time scale, which charac- 
terizes the exponential law, is Maxwell's time of relaxation, T^/j . In general it is defined by 

Tm = -. (7) 

cr 

In particular, in region 1, one has Tjv/i = ei/ci. In the limit of large time, the continuity equation implies 
V • j = 0, i.e. there are no sources or sinks. Also in this limit the electric potential satisfies Laplace's 
equation, ISV = 0, which follows from Gauss' law. Maxwell's relaxation time Tm^ in region 1 is very short, 
of the order of 10"^" s, thus the charge density tends very quickly to zero, and so does the extracellular 
potential. 



D. Behavior in region 2 

Region 2 is the near neighborhood of the membrane of a passive cell. The neuronal source creates the 
source field (or primary field) E source- As pointed out above (Fig. 01^]), due to the presence of a source and 
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the presence of free charges near the passive cell membrane, there will be a polarisation of the free charges, 
described by an electric induced field (or secondary field) Eind also denoted by Efree- Let us assume that 
the source field is "switched on" at time t = 0. Such time-dependence is described by a Heaviside step 
bmction H{t) (see Eq. jS^] in the Appendix). 

Esourceix.t) = Eo{x) H{t) . (8) 

The electric field present in region 2 results from the source field E source, the field due to free charges Efree 
(which create the induced field) and the field due to fixed localized charges (dipoles) of the membrane 
Emembr- Undcr thc hypothcsis that the passive cell membrane is a rigid structure with dipoles in fixed 
locations and undcr thc assumption that ion channels in that membrane remain closed, wc conclude that 
the electric field E„iembr does not vary in time. Gauss' law and the continuity equation now read 

Y7 P — Pfree 
V • £jfree 



£2 

V.j = -^. (9) 



Ohm's law reads 

j — ^2{E source -^/ree Emembr') , C-^^) 

which implies 

— ^^K:^^-^ ■{Esource+ Efree +E„,embr)=~^+ fix) for t > . (11) 

While pfree dcpcnds on position and time, the function /(x) denotes a time-independent term {Emembr is 
time-independent and E source is time- independent for i > 0) . The solution of Eq. [llj becomes 

Pfree{x, t) ^ a(f) 4^ b{x) cxpht/TA/J for < > . (12) 

Here, = e2/o'2 denotes the Maxwell time in region 2. The function a{x) represents the free charge 
density at equilibrium, that is a long time (t = oo) after the source has been switched on and the free 
charges have settled in region 1 in such a way that in region 1 no net electric field is left and hence no 
flow of electric current occurs. Thc function b{x) denotes the difference of the free charge density at the 
moment when switching on thc source and the free charge density at equilibrium. One should note that 
Pfreeix,t) is a continuous function in a; at < = 0. Poisson's equation implies that a linear relation holds 
between the free charge density and thc induced potential. This and Eq. yields 

V.ndix, t) = c{x) + d{x) exp[-i/rMj for i > . (13) 

In order to understand thc meaning of thc functions c{x) and d{x) , let us consider as example thc following 
source potential 

Vsource{S!,t) = (3H{t) + a[l - H{t)] , (14) 

where a and /? are taken as constant in space and time and H{t) is the Heaviside function defined in 
the Appendix. This source function makes a jump immediately after t = 0. Then thc induced potential 
becomes 

V.r^4x,t) = C7'(x) + [V.Mx,t^ +S) VZ"'ix)]oxp[~t/TM.] for t > , (15) 



where V^^^^^{x) denotes the induced potential at equilibrium and Vindix^t = +5) denotes the induced 
potential immediately after the source has been switched on. The resulting total potential is given by 

V{X, t) = Vsource{x, t = +5) + V^nd{x , t) for t > , (16) 

where Vsource {x, t = +5) denotes the source potential, Eq. |14j , immediately after the source has been 
switched on. 

Now let us consider more specifically a source with timc-dcpcndence given by the following function, 

V case(a) , , . 

Vsource(X,l) " | ^[^ _ ^(i)] case(6) . ^ ^ 

In case (a), the source potential is constant (= /3) in space and jumps in time from zero to one immediately 
after t = 0, meaning the source is switched on. In case (b), the source potential is constant (= a) in 
space and jumps in time from one to zero immediately after t = 0, meaning the source is switched off. By 
repetition of alternate switch-on's and switch-off's, one can introduce a temporal pattern with a certain 
frequency, as shown in Fig. |5J^] . Here we arc interested in the temporal behavior of the induced potential 
under such circumstances. In case (a) Eq. ^3] implies for the induced potential the following relation 

V^nd (.T, t) = V;;^' (a^) [1 - exp [-t/TMj] for t > . (18) 
By differentiation, one finds that the induced potential obeys the following differential equation, 

^Y^I^ + -^V^Ux, t) = ^VZ'\x) . (19) 

Similarly, in case (b) one obtains for the induced potential 

V,nd{.x, t) = V^ndix, t = o) cxp[-i/TA/,] for i > , (20) 
which obeys the following differential equation, 

dV,,nd{x,t) 1 



dt Tm, 

Eqs. |19I2H can be expressed as 



V,„d{x,t)^0. (21) 



dVind{x,t) , 1 1 rf^^rJU\ 

Clt 1 1 M2 

f(x)~ ( ^mr'(^) case(a) 

■'^^ \ case(6) ' 



(22) 



In order to solve this differential equation, one has to know the function f{x) given by V^^^^^'' (x) in case 
(a). How to obtain this function? We consider the case of an ideal dielectric medium in region 1. The 
function T^;^'/* (x) is related to the full potential and the source potential via 

= KT'' - KXK^) (23) 

The resulting total potential at equilibrium T^ot"'' on the surface of the membrane will be constant in time 
and in space (on the membrane surface), i.e. independent of position x, because otherwise there would be 
a flow of charges. The value of V^o^"' can be found by computing 

dS-e,Etr'=Q2f (24) 
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The integral is done over any surface S in the extracellular fluid, chosen such that it only includes the 
passive cell. Qlo™'^ is the total charge in the interior of such surface S. The physical value of Vjo^"'' must 
be such that the corresponding electric field Efot" yields QtoT — '^i^' Eq- , because the total charge 
of the passive cell before switching on the source was neutral. 

Note that the above reasoning is valid for an ideal dielectric medium, which is not the case for extracellular 
media. However, the small amplitude of the currents involved (~ 100 pA), the value of conductivity of 
extracellular fluid 3.3 S/m), the small dimension of most passive cells 10 fim diameter; ~ 2 nm of 
membrane thickness), and the high resistivity of membranes, imply a weak voltage drop on cell surfaces 
due to the current. Thus, the electrostatic induction is very close to an ideal dielectric. 



E. Source given by periodic step function: an equivalent RC electric circuit 

We considered above initial conditions where the source potential has been switched on at some time. 
This can be directly generalized such that the temporal behavior of the source potential is given by a 
periodic step function, with a period T5, as shown in Fig. The corresponding induced potential, 

presented in Fig. |1J3], shows a piecewise exponential increase followed by an exponential decrease. The 
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FIG. 2: Time-dependence of induced electric potential in response to an external source given by a periodic 
function. A. Source potential given by a periodic step function H{sin{u}t)), with period Ts = 2Tr/u> = 80ms. B. 
Induced electric potential for various values of Ts/Tm- In the case Ts/Tm « 1 (top), the induced potential 
fluctuates closely around the mean value of the source potential. In the case Ts/Tm >> 1 (bottom) the induced 
potential relaxes and fluctuates between Vmin = mV and Vmax ~ 100 mV of the source potential. As a result, 
the amplitude of oscillation becomes more attenuated with increased frequency / = l/Ts- 
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figure shows the response to a source of periodic step function, for different values of Tg measured in 
units of Tm- One observes the following behavior. At the top of the figure the period Ts is smallest, i.e. 
which corresponds to a rapid oscillation about its time average. At the bottom of the figure the period 
Ts is largest, corresponding to a low-frequency oscillation. It is important to note that the amplitude of 
oscillation (i.e. the difference between its maximal and minimal value) is much smaller for small Ts than 
for large Ts- This yields an attenuation effect of high frequencies in the induced potential. 

Such behavior of the induced potential is well known from and mathematically equivalent to that of an 
RC electric circuit, with a resistor R and a capacitor C (see Fig. |3]). The equation of motion in such an 




wv 



source 



A 



V 



ind 



V 



FIG. 3: Model of electric circuit with capacitor C, resistor R. Vsource denotes the potential of the source and 
denotes the induced potential at the capacity C. 

RC circuit relating the induced potential V^^j^ to the capacity C to the source potential Vsource is given by 



dt RC ™ ^ source\i-) 



(25) 



This equation is mathematically equivalent to Eq. 



if we identify 



equil 



ind 



G 

€ ^ 

Tm 



x)H(t) 

^R, 
C , 

^ RC . 



.it), 



(26) 



Thus, for each point x near the membrane of a passive cell (region 2 in Fig. 0A.]), we can set up an 
equivalent electric circuit, which is equivalent in the sense that it gives a differential equation with the 
same mathematical solution as for the induced potential of the original system of a source and a passive 
cell. 



F. Induced potential for time dependent source 

The analogue of an RC electric circuit can be used to determine the characteristic properties of the 
source/passive-cell system, like the cut-off frequency. For that purpose let us consider the source potential 
to be given by a periodic function of sinusoidal type switching on at time t = 0: 



source 



{x,t) = Vo{x) exp(iwt) H{t) 



(27) 
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The asymptotic (large t) solution of the induced potential is given by 

Vind{x,t) ^t-.oo — r~- — Vo{x) exp{iujt) . (28) 
(7 + iLue 

Fig. [33] shows an example of the induced potential, where the source has a time dependence given by 
an oscillating step function (Fig. [2K])- We use the notion of the transfer function, which denotes the ratio 
of induced potential corresponding to frequency / of the source over the induced potential at frequency 
/ = of the source. 

Here, V^^{x, t) denotes the solution of the differential equation [SJ with a source term having an asymptotic 
time-dependence given by exp(ia;<). There is a variety of methods to analyze the time-dependence of 
induced potentials for a given time-dependent source. The method of Fourier analysis, where the time 
dependence of a function is decomposed in terms of function exp(ia;t), has been proven very useful in the 
time-analysis of experimental data (e.g. the power spectrum). An alternative method is the Heaviside 
method, where the time dependence of a function is decomposed in terms of oscillating step functions 
H{sm{ujt)). In the Appendix we discuss the asymptotic behavior of the induced potential, where the 
source term is given in terms of Heaviside step functions and also in terms of Fourier components. Then 
the transfer function for a source given by a Fourier component exp(iwt) is given by 

Ftm{^) = ■ (30) 

Ftm is a complex function, with its modulus given by 

I^tm(-)| ^ ] ■ (31) 

For frequency zero one has |Fta/| = 1. The cut-off frequency fc is defined such that falls off to the 

value |Fta/| = l/\/2. Thus we find 

/c = g = ^ = (2.r..)- . (32) 

The modulus of the transfer function vs. frequency is shown in Fig. 0], for the case of parameters a and 
e of region 2, as given in sect. ljll A|l . This corresponds to a cut-off frequency fc = 100 Hz. Such behavior 
represents a frequency filter. It shows that we can compute the induced potential as a function of frequency 
based on the analogy with an equivalent electric circuit with a resistor R and capacity C arranged in serial 
order. A series of detailed numerical simulations, modeling frequency filtering of extra-cellular neural tissue 
in the neighborhood of passive cells in terms of the electric circuit model, Eq. j^S], using physiologically 
realistic parameters of e and a are shown in Fig. [3]. They present the signal coming from a periodic step 
function with a period Ts = 10 ms, i.e. a frequency of 100 Hz. Going from Fig. [5j^] to Fig. [Sp] T5 is 
kept fixed but Tm increases, i.e. Ts/Tm decreases. Fig. corresponds to the case where the relaxation 
time T^v/ is 10 times smaller than the period of the signal Ts- In Fig. |SJ3] Ta/ and Ts are equal, and in 
Fig. [Sp] Tm is 10 times larger than Ts- We observe that the shape of the signal after filtering is more 
or less intact for Ts/Tm = 10 and becomes gradually more deformed when going over to Ts/Tm = 0.1. 
As the Fourier transformation after filtering shows, the higher frequency components get gradually more 
suppressed. This represents a low frequency band pass. 

To close this section, we would like to remark that the treatment above is general and not restricted 
to step functions. The Heaviside functions were used as a tool to distinguish transient and asymptotic 
behaviors, as well as to calculate the transfer function. For an arbitrarily complex time-dependent source, 
the convolution integral with this transfer function gives the filtered signal, so this approach can be applied 
to any physical signal. 
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FIG. 4: Modulus of transfer function vs. frequency, using biologically realistic parameters e and a. 

III. MODELING FREQUENCY-DEPENDENCE OF MULTIPLE NEURONS AND PASSIVE 

CELLS 

In Section [H] , we laid out the theoretical framework of the model and showed that this model ~ with 
respect to time and frequency-dependence but neglecting spatial dependence - is equivalent to a RC-circuit 
model. As an example, we have presented numerical simulations treating the case of a single source and 
a single passive cell. In this section, we would like to to generalize this approach to take into account 
space- and time-dependence in LFPs. In this and the following section, we present numerical simulations 
considering more complex three-dimensional arrangements of neuronal sources and passive cells. 

A. Computation of field potentials for multiple neurons and passive cells 

In a first set of numerical simulations, we considered three-dimensional arrangements of sources and 
passive cells. To solve the equations, we used a discretization method similar to finite elements methods (see 
below). The neuron source was represented by a cubic shape of length Lneuron, and likewise, we consider 
the passive cells were represented by cubes of length Lpass- Neuron source and passive cells were located 
at some distance d from each other. Although being far from realistic, this geometry facilitates numerical 
simulations. Since the aim of this work is not to make quantitative predictions, such an arrangement will 
help explaining qualitatively the frequency attenuation properties of field potentials in three-dimensional 
arrangements of cells, as we will show below. 

Let us first outline the method of numerical computation of the potential. From a mathematical point of 
view, the solution of a partial differential equation requires to specify suitable boundary conditions. In the 
context of Maxwell's equations, Dirichlct's boundary conditions (potential given on a closed surface) or von 
Neumann's boundary conditions (electric field given on a closed surface) arc known to guarantee a unique 
solution. Here, we suggest to use boundary condition of Dirichlet's type, that is we assume the potential 
is given on the surface of all cells. The problem is, however, that the value of the physical potential on the 
surface of those cells is not known a priori, and, moreover, the functional values of the potential can not 
be chosen freely. The physical solution is determined by the principle that all free charges are located such 
that the electromagnetic field energy becomes minimal (Thomson's theorem). This imposes constraints on 
the boundary conditions. 

How can we figure out those a priori unknown boundary conditions? We use the physical principle 
that those physical boundary conditions arc equivalent to the condition that Gauss' integral summing the 
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FIG. 5: Effect of filtering of a periodic electrical potential representing a neural source. A. The original signal 
(top row, left panel). The frequency spectrum of such signal is obtained by Fourier transformation (top row, right 
panel). Then the signal is sent through a filter (middle row, right panel). After an inverse Fourier transformation 
the filtered signal as function of time is obtained (middle row, left panel). The period of the original signal is 10 
ms, and the relaxation time of the filter is 10 ms. B. Same signal as in A, filtered using a relaxation time of 1 ms. 
C. Same signal as in A, filtered using a relaxation time of 100 ms. 



flux around a charge source becomes invariant when the integration surface is changed. A computational 
strategy to find those boundary conditions is to use a variational principle, starting from some initial guess 
of boundary conditions and then to adjust iteratively the boundary conditions until eventually the total 
electric energy attains its minimum. 

We used a variational principle via the following iterative scheme. In a first step, we compute the 
potential due to the source in absence of any passive cell in the asymptotic regime (a long time after the 
source has been switched on), for the case where the source has no frequency dependence {lo — 0). We 
assume that the potential obeys the boundary condition on the surface of the neuron source 

V{x) = = const for x on surface of source . (33) 

The potential satisfies Laplace's equation 

AV{x) = . (34) 

The numerical solution has been carried out using the standard method of relaxation. 

In a second step, we place a passive cell at some distance d from the neuron source. We want to calculate 
the induced potential due to the presence of the passive cell. We proceed in the following way. Let us 
consider for a moment that the passive cell is small compared to the size of the neuron, Lpass << ^neuron- 
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Then the induced potential on the surface of the passive cell would be almost identical to Vp^'!}^'^^ , the 
potential created by the source alone, evaluated at the center of the passive cell. Then imposing the 
boundary condition that the potential on the surface of the passive cell takes the value Vp^^g*^^ would be 
very close to the exact solution. Now wc consider the passive cell no longer as tiny. Then imposing the 
boundary condition boundary on the surface of the passive cell 

V{x) = Vp^'l^^'^^ = const X on surface of passive cell (35) 

becomes an approximation. 

In a third step, taking into account that the potential obeys the boundary conditions on the neuron 
source and on the passive cell, Eqs. j33l35| . wc solve again Laplace's equation, Eq. [33]. Now we test if the 
obtained solution V{x, t) is the correct solution. Using Gauss' theorem, one has 

/ ds-E^ I ds-E . (36) 

Here the integral denotes an integral over a closed surface, first englobing the neuron source and second 
englobing the neuron source plus the passive cell. Because the passive cell has a total charge Qpass = 0, 
Eq. should be satisfied by the physical solution Vphys{x,t). As long as V{x,t) differs from Vphys{x,t), 
one has to adjust the potential on the surface of the passive cell, solve again Laplace's equation and verify 
the charge balance Eq. |3t)j . This a an iterative process, which turned out to converge quite fast towards 
the physical solution. 

This method can be applied as well to treat multiple neuron sources and multiple passive cells. We would 
like to point out that the extracellular neural tissue is composed of much more passive cells than neurons. 
On average the number of glial cells per unit volume is higher than the number of neurons per unit volume 
by about a factor of 10. The physical consequence of this property is that the induced potential due to the 
presence of passive cells becomes more important than the source potential. 

Finally, to investigate the effect of time-dependent sources and multiple passive cells, we take advantage 
of the following equivalence: the induced potential on a cell with a relaxation time Ta/ subject to a time 
dependent source is equivalent to the induced potential in a cell with Tm = subject to a "filtered" source 
given by the convolution of the source with the transfer function [Ftm in Eq. [21]). This corresponds to the 
following calculation steps, for each passive cell: (1) evaluate the source S{x, uj) provided by the neighboring 
system of sources and other passive cells; (2) evaluate the "filtered" source S*[x, u) = S[x, uj)Ftm{i^) (from 
Eq. 1201); (3) calculate the induced potential on the passive cell using a similar procedure as outlined above 
for a time- independent source (w = 0). The whole procedure is repeated iterativcly until all sources and 
induced potentials converge towards the physical solution. Numerical results using this iterative procedure 
are presented in the next section. 

It is important to note that in the range of frequencies considered here, the wavelengths are too large 
for significant wave propagation phenomena. Accordingly, we only consider sources where space and time 
dependences factorize. Consequently, the transfer function does not have any space dependence, as shown 
in the calculation in the Appendix, Eqs. |51l66j . In addition, the application of the transfer function holds 
for spherical sources as well as for cubic sources. 

B. Numerical results 

We now present results of numerical solutions of Laplace's equation for the scenario in the presence 
of a source and several passive cells. We have taken as parameters of conductivity and permittivity 
those corresponding to the intracellular space at the membrane of the passive cell, given in sect. |II A| 
(conductivity is non-zero). It corresponds to a cut-off frequency of fc = 100 Hz. The corresponding 
Maxwell time and transfer function arc obtained by Eqs. |32l30j . 
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The goal is to show that the induced potential gives an important contribution to the total, i.e. physical 
potential. The results are presented in Figs. [nHH]. For the case of a single source and a single passive cell 
the source potential, the induced potential and the total potential are shown in Fig. HK-C], respectively. 
Comparing the source potential (B) with the induced potential (C), one observes that the latter gives 
a substantial contribution. In Fig. 0\-D] we show the induced potential multiplied with the norm of 
the transfer function. We observe that the frequency-dependent induced potential decreases when the 
frequency of the source increases from to 400 Hz. Figs. |SH5] show the corresponding results in the case 
of a single source and four passive cells. We observe qualitatively the same results as in Fig. [S]. These 
results demonstrate clearly the phenomenon of attenuation of high frequencies in the induced potential. It 
should be mentioned that the computations were carried out in 3-D, but are depicted as projection in 2-D. 
Qualitatively the same results are obtained in the case of two sources and a single passive cell (not shown) . 
Again we found for the time-independent source that the induced potential is non-negligible compared to 
the source potential and the frequency-dependent induced potential rapidly decreases when the frequency 
goes from to 400 Hz. 




FIG. 6: Extracellular potential generated by a time-independent source situated nearby a passive cell. Both the 
source (s) and the passive cell (a) have a cubic shape; the passive cell was twice smaller than the source. A. 
Total extracellular electric potential (Vsource + Vi„d) in a horizontal plane cut at the base of both cells. B. Electric 
potential resulting from the source only (Vaource). C. Induced electric potential (Vind)- AH potentials were obtained 
by solving Laplace's equation with a static source (/ = 0), and are represented in units of percent of the source 
potential. The distance is represented in units of 0.1 fim. 
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FIG. 7; Same as Fig. but showing the induced potential Vind for a time-dependent source (frequencies / = 
0, 100, 200,400 Hz in A,B,C,D, respectively). The induced potential tends to zero when / increases. 

IV. ATTENUATION AS A FUNCTION OF DISTANCE 

The simulations shown in the preceding section illustrated the fact that the induced potential vanishes 
for very high frequencies of the source field, a fact that can also be deduced from Eq. [22] . In other words, 
for very high frequencies (>> T^/), the extracellular field will be equal to the source field, since the induced 
field will vanish. The space dependence is easy to deduce in such a case, and the extracellular potential 
attenuates with distance according to a 1/r law, as if the source was surrounded only by conducting fluid. 

However, for very low frequencies (<< T^^), the space dependence of the extracellular potential will 
be a complex function depending on both the 1/r attenuation of the source field, and the contribution 
from the induced field. Such a space dependence is not easy to deduce, since it depends on the spatial 
arrangement of fluids and membranes around the source. In this section, we attempt to derive such a low- 
frequency space dependence for a more realistic system of densely packed cells (illustrated in Fig. [TUK]'). 
To constrain the behavior to low frequencies, we only consider the zero-frequency limit by using a constant 
source field. We proceeded in two steps. First, wc calculate the electric potential at the surface of a passive 
cell (Section |I V Aj ) . Second, wc calculate the spatial profile of LFPs in a system of densely packed spheres 
of identical shape (Section |IV B| ) . 
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FIG. 8: Extracellular potential generated by a time-independent source surrounded by six passive cells. Same 
description as in Fig. but for a system of a single static source (s) and 6 passive cells. The figure shows the 
potential in a horizontal plane in which 4 passive cells are visible (a). A: Total extracellular potential, B: Source 
electric potential, C: Induced electric potential Vind- 



A. Electric potential at the surface of passive membranes at equilibrium 

Let us assume a spherical passive cell embedded in a perfect dielectric medium, and exposed to a constant 
electric field. At equilibrium, we have seen above that the effect of the electric field is to polarize the charge 
distribution at the surface of the cell, such as to create a secondary electric field (see Fig. ^3]), but the 
induced electric field is zero inside the cell. In this case, the conservation of charges on the surface implies: 



psurf dS = , (37) 

Surf 

where psurf is the charge density on the surface of the cell. The resulting electric potential is given by: 

Vtot{x,y,z) =Vsource{x,y,z) + PSMr/ ^g, ^^g^ 

J J Surf 47re r 
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FIG. 9: Same as Fig. but showing the induced potential Vi„d for a time dependent source at different oscillation 
frequencies (/ = 0, 100, 200, 400 Hz in A,B,C,D, respectively). 

where Vsource is the electric potential due to the source field, Vtot is the total resulting electric potential 
due to the source field and the induced field, and r is the distance from point (x,?/, z) to the center of the 
cell. Because at the center of the cell, (a, 6, c), we necessarily have r = R (where R is the cell's radius), the 
value of the resulting electric potential at the center is given by: 

Vtot{a, b, c) = Vsource{a, b,c)+ / / dS = Vsource{a, b, c) . (39) 

J J Surf 47re K 

Thus, the electric potential at the surface of a spheric passive cell at equilibrium equals the potential due 
to the primary field at the center of the cell. In other words, the effect of the secondary field in this case 
perfectly compensates the distance dependence of the primary field, such as the surface of the cell becomes 
isopotential, as discussed above. 

B. Attenuation of electric potential in a system of packed spheres 

Keeping the assumption of a constant electric field, we now calculate how the extracellular potential 
varies as a function of distance in a simplified geometry. We consider a system of packed spheres as 
indicated in Fig. [TUR]. From the solution of Laplace's equation, the extracellular potential at a distance r 
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FIG. 10: Extracellular potentials as a function of distance for a system of densely packed spherical cells. A. Scheme 
of the arrangement of successive layers of identical passive cells packed around a source (in gray) . The potential of 
the source is indicated by Vq. Vi, V2 ... Ki indicate the potential at the surface of passive cells in layers 1,2, ...n, 
respectively. The dotted lines indicate isopotential surfaces, which are concentric spheres centered around the 
source, and which are indicated here by Si, S2 ... 5„. B. Extracellular potential as a function of distance (in units 
of cell radius), comparing two cases: with induction (solid line, corresponding to the arrangement schematized in 
A), and without induction (dashed line; source surrounded by conductive fluid only). Both cases predict a different 
scaling of the electric potential with distance (see text for details). 



from the source center is given by: 

V{r) = - , (40) 
r 

where is a constant, which is evaluated from the potential at the surface of the source (5*0 in Fig. |l()l\]): 

k 



ViR) - Vo 



R ' 



(41) 



where R is the radius of the source. Thus, the potential due to the source field at a given point r in 
extracellular space is given by: 



V{r) 



RVo 



(42) 



Considering the arrangement of Fig. jlOl^]. if all cells of a given layer are equidistant from the source, 
their surface will be at the same potential (sec Section |IV A| ). which we approximate as a series of 
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isopotential concentric surfaces {Si, S2 ■■■ in Fig. jlUI^]'). A given layer (n) of isopotential cells is 
therefore approximated by a new spherical source of radius r„, which will polarize cells in the following 
layer {n + 1). According to such a scheme, the potential in layer n + 1 is given by: 

Vn+l = ^ , (43) 

where dn+i is the distance from the center of cells in layer {n + 1) to the center of the source. According 
to the scheme of Fig. |10l\]. we have r„ = {2n + 1)R and dn = 2nR. Thus, we can write the following 
recurrence relation: 

K+i = ^ . (44) 



Consequently: 



which can be written, for large n: 



^"+1 = (IT ^) ^0 , (45) 



(2n + l)! 2(2n)! 
^" - 22"(n + l)!n! ^° ¥^ ^° ' 



Using Stirling's approximation, n\ ~ (n/e)"'\/27m for large n, leads to: 



y„ ~ ^ . (47) 



Thus, in a system of densely packed spherical cells, the extracellular potential falls-off like 1 / ^ (Fig. jlOb]. 
continuous line). 

In contrast, in the absence of passive cells in extracellular space, the electric potential is given by the 
source field only (Eq. |42j ) . which, using the same distance notations as above, is given by: 

Such l/r behavior is illustrated in Fig. |lUb] (dashed line). Note that other theories predict a steeper 
decay. For instance the Debye-Huckel theory of ionic solutions predicts a fall-off as exp(— fcr)/r. 

Thus, for this particular configuration, there is an important difference in the attenuation of extracel- 
lular potential with distance. The extracellular potential in a system of densely packed spheres falls-off 
approximately like in contrast to a l/r behavior in a homogeneous extracellular fluid. Note that a 

behavior can also be found in a system in which the source is defined as a current. For a constant 
current source /g, with variations of conductivity following a spherical symmetry around the source, the 
extracellular potential is given by [Tsj : 

V{r) = / / ^^dr, (49) 
47r r^a{r) ^ 

where CT(r) is the radial profile of conductivity around the source. Assuming that a{r) ~ ao/y^, gives 

V{r) - = J^Vo. (50) 

Consequently, the l/y/r behavior found above is functionally equivalent to a medium with conductivity 
varying like 1 / \/r . This "effective conductivity" is similar to that introduced in a previous study . 
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V. DISCUSSION 

From theoretical considerations and numerical simulations we have obtained the following main results: 
(i) We explored the assumption that neuronal current sources polarize neighboring cells. This polarization 
produces an induced electric field, which adds to the electric field directly produced by the sources. This 
induced field is non-negligible for biologically realistic parameters, (ii) This induced electric field has 
strong frequency dependent properties. This system is equivalent to an equivalent RC circuit and always 
has low-pass filtering properties, (iii) The cut-off frequency of this low-pass filter is determined by Maxwell's 
relaxation time of the membrane surfaces surrounding neuronal sources, (iv) Consequently, the attenuation 
of high frequencies will not be influenced by this induced field, and will be the same as if neurons were 
embedded in a homogeneous extracellular medium. On the other hand, the attenuation of low frequencies 
will depend on the induced field, and will therefore depend on the geometry of fiuids and membranes in 
the extracellular space. 

In previous work, we showed that inhomogeneities of conductivity or permittivity in extracellular space 
could give rise to strong frequency filtering properties 0. However, low-pass or high-pass filters could be 
obtained depending on the profiles of conductivity used in this model. In the present model, a low-pass 
filter is predicted, consistent with experiments (which never evidenced a high-pass filter). However, the fact 
that frequency filtering properties are due to the alternance of high-conductive fiuids and low-conductive 
membranes in the present model is compatible with the inhomogcneity of conductivity postulated in the 
previous model. One main difference with the previous model is that here, we explicitly considered a 
non-zero charge density on neighboring membranes. 

Another similarity with the previous model concerns the mechanism of attenuation of LFPs with distance. 
The primary or source field experiences a steep attenuation for all frequencies, and the electric potential 
will attenuate with distance r following a l/r law, similar to the attenuation in a homogeneous conducting 
fiuid 0]. On the other hand, the induced field will be subject to a different law of attenuation with distance, 
which will depend on the spatial arrangement of pass ive cells around the source. In the previous model 
it depended on the particular conductivity profile [ig. In general, for densely packed cellular membranes 
around the source, one will have a law of attenuation which will be of l/r", where a < 1 may depend on the 
frequency w. For example, we found a ~ 0.5 for spherical cells in Fig. jlO| ) for low w. Thus, both models 
predict the following scenario: high frequencies follow an attenuation in l/r, similar to a homogeneous 
extracellular fluid, but low frequencies follow a slower attenuation profile, because these frequencies are 
"transported" by the induced field. 

It should be noted that this model represents a strong approximation of the actual complexity of the 
mechanisms involved in LFP generation. A first approximation was to consider passive cells as spheres, 
neglecting their morphological complexity. Another assumption was that the electric field results only 
from neuronal current sources, while other cell types, such as glial cells, also contains ion channels and 
may influence LFPs. A third approximation was to neglect the effect of variations of extracellular ionic 
concentrations (like potassium buffering by glial cells), which may also influence LFP activity, especially at 
low frequencies (in addition to local variations of conductivity). It is difficult to estimate the consequences 
of these different approximations, but the polarization effects described here should occur in more complex 
situations, so it is likely that the present considerations apply to more complex geometries and current 
somces. This should be tested by numerical simulations, for example by constraining the model with 3D 
reconstructions of extracellular space from electron microscopic measurements. 

The scenario outlined above of a "transport" of low frequencies can be tested experimentally in different 
ways. First, measuring the decay of specific frequency components of LFPs with distance, and in partic- 
ular how they differ from the l/r law, should yield direct information on how much extracellular space 
deviates from a homogeneous conducting fiuid. The attenuation profile should also depend on the spatial 
arrangement of successive fluids and membranes in the extracellular fluid. It is conceivable to inject sinu- 
soidal currents (of low amplitude to allow induction effects) of different frequencies using a microelectrode 
and measure the LFP produced at increasing distances, to obtain the law of attenuation as a function of 
frequency (i.e., a{Lj)). The comparison with model predictions should tell us to what extent the model 
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predicts the correct attenuation. 

The equivalent RC circuit investigated above also deserves some comment. It is well known that an 
electric circuit consisting of a resistor and a capacitor, which stores electric energy, introduces a phase 
difference between current and potential, and has characteristic low-pass filtering properties. The fact that 
such a simple RC circuit may be used as a model of frequency filtering of cerebral tissue is, however, more 
surprising. In fact, if we construct an arrangement of densely packed passive cells around a current source 
(such as in Fig. |10t\]'l. the system will be equivalent to a series of RC circuits, each representing one 
"layer" of cells. In this context, it may be interesting to relate the type of spatial arrangement with the 
type of equivalent circuit (s) obtained, for example by considering different cases such as cells or random 
diameters, various shapes, or even by using 3-D morphological data from real brains obtained by serial 
reconstructions from electron microscopy. This type of investigation constitutes a possible extension of 
the present work. In addition, we have considered a source being given by a single neuron. In reality, the 
source of LFP activity is a system of many neurons with complicated phase relationships between their 
activities. An extension of our model in this direction would also be a most important step to take. 

A second model prediction that should be testable experimentally concerns the predicted cut-off fre- 
quency This information can be related to experimentally recorded LFPs which could be analyzed using 
Fourier analysis to yield information about the cut-off frequency. However, such an analysis should be done 
by comparing different network states, to make sure that the cutoff frequency is not dependent on activity 
but rather depends on structural parameters. If correctly done, this analysis should provide an estimate of 
the Maxwell's relaxation time for the surface of membranes, or equivalently the tangential conductivity of 
membrane, for which there exists presently no experimental measurement. Work in this direction is under 
way. 

Finally, we would like to conclude by suggesting that these theoretical predictions in conjunction with 
measurements of frequency-dependent field potentials may lead to new ways for detecting anomalies, for 
example due to degenerative diseases like Alzheimer or brain tumors. If such diseases cause structural 
changes in brain tissue, e.g. creation of vacuoles, this will manifest itself by changes of electric proper- 
ties like the attenuation of low-frequency components with distance. These changes should be visible in 
measurements of EEC, ECoG or LFP signals from multielectrodcs in the area of affected brain tissue. 
Detecting such changes supposes a prior understanding of the spectral filtering properties of these electric 
signals, which is one of the main motivation of the present work. 



Here we discuss the Heaviside method as well as the Fourier transformation for the purpose to determine 
the asymptotic behavior of the solution of a class of differential equations. We want to determine the 
asymptotic behavior of the following differential equation. 



VI. APPENDIX 



+ -V(x,t) = -B(x)H(t)cxp(iLot) . 

at t e 

The Heaviside function H{t) or 0-function 6{t) is defined as 




(51) 




(52) 



Let us start out by considering the following differential equation. 



(53) 



One verifies that the solution for t > Q has the form 



V{x,t) = B{x)[l 



exp( 1)] . 



(54) 
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For t — !■ oo, one finds 

V{x,t) ^t^oo B{x) . (55) 

That is B{x) is the boundary condition at t = oo for the solution of Eq. j^. On the other hand, if we 
consider the differential equation 

^^("'*)+^y(x,t)=0, (56) 



dt e 



and impose the boundary condition 



one obtains the solution 



V{x,t^O) ^ B{x) , (57) 

V{x, t) = B{x) exp(- ^i) . (58) 

Now let us consider as inhomogcneous term an oscillating step function, 

i7(sin(wt)) . (59) 

This function hops between zero and 1 with a frequency ta thus giving a oscillating step function. A 
function of this kind (however shifted and stretched) is shown in Fig. Now let us consider Eq. |55] 

with solution given by Eq. j54| . The inhomogcneous term of the differential equation (r.h.s.) is constant for 
positive time. The solution increases with time to approach an asymptotic value. Next consider Eq. j56| 
with solution given by Eq. [SS]. The inhomogeneous term of the differential equation (r.h.s.) is zero for 
positive time. The solution decreases with time to approach the asymptotic value zero. Thus, when we 
consider now the differential equation 

^Xi^ + 'Lv[x, t) = -B{x)H{t)H{sin{Lut)) , (60) 
at e t 

we see that its inhomogeneous term is an oscillating step function. That means its inhomogeneous term 
alternates between the inhomogeneous term of Eq. and that of Eq. jSS] . Consequently, the solution of 
Eq. inO] alternates between that given by Eq. jSl] and that by Eq. [SS], i.e. its increases and then decreases 
exponentially. Such a solution is shown in Fig. |33]. 

Now let us consider the differential equation |5H . We want to compute its asymptotic behavior for large 
positive time. We consider the homogeneous equation 

du{x,t,to) , cTt./ , , N n fd.s 
— \- -V{x,t,to) = . (61) 



It has a solution 



u{x,t,tQ) = B{x) cxp{iu}tQ 1) , (62) 



which satisfies 

u{x, t = 0, to) = B{x) exp(iwto) • (63) 

Based on the solution u{x,t,tQ), Eq. of the homogeneous Eq. one can express the solution of the 
inhomogeneous Eq. |5H by 



V{x, ^) ~ J B{x) exp (j,ujtQ — — (t — to)j 

= exp(ia;t) / dr _B(a:) exp(— iLjr r) . (64) 

Jo ^ 
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Its asymptotic behavior is given by 

V{x,t) ~t-»oo cxp{iujt) / dr B{x) cxp{—iujT t) 

Jq £ 



Then from the definition of the transfer function, Eq. |29j . one obtains 



^—^ B{x) cxp{iujt) . (65) 



Ftm{uj) = T— ^ . (66) 
1 + luje/a 
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